Degree distributions in mesoscopic and macroscopic functional brain networks 
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We investigated the degree distribution of brain networks extracted from functional magnetic res- 
onance imaging of the human brain. In particular, the distributions are compared between macro- 
scopic brain networks using region-based nodes and mesoscopic brain networks using voxel-based 
nodes. We found that the distribution from these networks follow the same family of distributions 
and represent a continuum of exponentially truncated power law distributions. 

PACS numbers: 87.19.If, 02.10.Ox, 87.18.Sn, 89.75.Da 



Small-world networks are a class of networks char- 
acterized by highly interconnected neighborhoods and 
efficient long-distance connections, connecting any two 
nodes in a network with just a few intermediary con- 
nections fl*]. Since the introduction of small- world net- 
works, these small-world properties have been observed 
in many social, technological, and biological networks [2[. 
The network organization of the human brain has also 
been demonstrated as a small- world network [3, 4]; small- 
world properties have been verified in both anatomical 
B B a 11 and functional [1, [TO, lU, E, El, H brain 
networks. A small- world structure is advantageous for 
brain networks since it can support both localized pro- 
cesses specific to different brain regions as well as dis- 
tributed processes encompassing multiple brain regions 
at once. Some studies have also reported that functional 
brain networks are scale-free networks net- 
works characterized by mega-hubs with extremely large 
node degrees and by the degree distribution following a 
power law distribution [15]. 

Interestingly, scale-free properties have been observed 
in brain networks using voxels (3 p p ixels in 3D brain 
images) as network nodes [ll|, [l2|, [l3[, but not in brain 
networks modeled using anatomical regions as nodes 
[5|, 16|, 18|, [T^. The distribution from a region-based 
network is not truly scale-free but follows an exponen- 
tially truncated power law distribution [H, 0, 0]- This 
discrepancy is rather surprising since both types of net- 
works describe the same biological system, the human 
brain. Although the reason for the difference is unclear, 
it may be due to differences in data processing steps in 
constructing the brain network [4] . Another possible rea- 
son for the discrepancy is the difference in the scale of 
these networks; while voxel-based networks represent a 
finer mesoscopic organization of the brain, region-based 
networks represent a coarser macroscopic organization of 
the brain. Thus, to examine the effects of data processing 
and the network scale on the node degree distribution, we 
constructed region-based and voxel-based networks based 
on the same functional MRI (fMRI) data. We compared 
the degree distribution from both types of networks de- 
scribing the organization of the same set of brains. 

The data set for this study consisted of fMRI exper- 
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FIG. 1: A schematic of data processing steps to generate 
a region-based network and a voxel-based network from an 
individual fMRI data set. Voxel time courses are extracted 
from a series of 3D fMRI images. In a region-based network, 
extracted time courses are averaged in different anatomical re- 
gions to produce a node time course in each region, whereas 
in a voxel-based network each voxel time course is treated as a 
distinct node time course. A correlation matrix is generated 
between node time courses, and then thresholded by a cor- 
relation threshold Rc to produce a binary adjacency matrix 
representing a network. 



iment data from 5 subjects. For each subject, a series 
of MRI images measuring neurological activities were ac- 
quired every 2.5 sec. for 5 minutes, resulting in 120 im- 
ages, each with 46 x 55 x 37 voxels of size 4 x 4 x 5mm. FIG 
[1] displays the schematic of the data processing steps. To 
construct a voxel-based network, each of ~20,000 brain 
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voxels was treated as a node. For a region-based net- 
work, the voxel time courses were averaged for voxels 
within each anatomical region, resulting in time courses 
corresponding to approximately 100 nodes of anatomi- 
cal regions. In both networks, the correlation coefficient 
was calculated between the time courses of two nodes, 
producing a correlation matrix. If the correlation coef- 
ficient exceeded a correlation threshold Rc^ these nodes 
were considered to be functionally connected. By ap- 
plying a threshold Rc to the correlation matrix, a bi- 
nary matrix was formed, known as an adjacency ma- 
trix, with 1 indicating existence of an edge connecting 
two nodes and otherwise. To facilitate a comparison 
of degree distributions across subjects, Rc was adjusted 
for each subject so that the average node degree for the 
network was similar across subjects. The average de- 
gree was 30.3± 1.3(SD) for the voxel-based networks and 
4.44 ± 0.19 for the region-based networks. 

FIG [2] show the degree distributions resulting from the 
voxel-based and region-based networks. FIG [2](a) shows 
the degree distributions of the voxel-based networks plot- 
ted on a log-log scale. For all the subjects, the degree 
distribution P{k) seems to follow a straight line as seen 
in other voxel-based networks , indicative of a 

power law distribution P{k) oc with 7 1.4. How- 
ever, the tail of the distributions exhibit increased uncer- 
tainty, giving an appearance of a fuzzy tail. To under- 
stand the distribution profile better [16], we calculated 
the cumulative distribution F{k) = T,k'<kP{k') and plot- 
ted the complementary cumulative distribution 1 — F{k) 
on a log-log scale in FIG[2]^b). If the distribution were 
truly a power law distribution, then the plot of 1 — F{k) 
would also follow a straight line [16]. Instead, the distri- 
butions decay faster than a power law distribution, and 
follow an exponentially truncated power law distribution 
P{k) (X kf^ exp{-k/kc) with p = -0.61 and kc = 81.3 
(dashed curve, FIG[2](b)). 

FIG[2](c) shows the log-log plots of 1 — F{k) from the 
region-based networks. The distributions are curved, 
showing an accelerated decay for higher k. These 
distributions follow a gamma distribution P{k) ex 
/c*^"^ exp(— /c/Zcc), which is another parameterization of 
an exponentially truncated power law distribution with 
the exponent a — 1 instead of /3. As a reference, FIG[2fc) 
also shows the best-fit gamma distribution, with a = 2.87 
and kc = 1.60 (dashed curve). These parameters are 
within a similar range as other region-based functional 
and anatomical brain networks [sl, @, B S [l3] • 

The results from both types of networks demonstrate 
that the degree distribution follows an exponentially 
truncated power law distribution. This is not surpris- 
ing since both networks are modeling functional connec- 
tivities of the same brains, although scales are differ- 
ent. As a comparison, we also fitted other distributions, 
namely a power law distribution P{k) cx k~^ and an ex- 
ponential distribution P{k) oc exp(— /c//cc), but neither 
approximated the distribution better than an exponen- 
tially truncated power law distribution. It is interesting 
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FIG. 2: Node degree distributions from the voxel-based and 
region-based networks: the probability distributions P{k) 
from the voxel-based networks (a), the complimentary cumu- 
lative distributions 1 — F{k) from the voxel-based networks 
(b), and the complementary cumulative distributions from 
the region-based networks (c). The best-fit curves of exponen- 
tially truncated power law distribution P{k) oc exp(— /c/A^c) 
(dashed curves) are also shown, with p = —0.61 and kc = 81.3 
for the voxel-based networks (b) and /3 = 1.87 and kc = 1.60 
for the region-based networks (c). 

that the degree distributions P{k) from voxel-based net- 
works appear to follow a power law distribution as pre- 
viously reported by other studies. However, the plots of 
the cumulative distributions 1 — F{k) indicate otherwise. 
Therefore the scale-free properties seen in other voxel- 
based networks [ll|, [H, [l^ may simply be an artifact of 
data display. Since there are very few nodes with high /c, 
the empirical estimation of P{k) based on a histogram is 
very unstable as noted by Keller fidl\ , and a straight line 
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TABLE I: The estimated exponent f3 and cut-off parameter 
kc for the exponentially truncated power law distributions for 
the networks with various scales, along with the number of 
nodes N and the average node degree (k). 



Network 


N 


(k) 




kc 


Voxel-based network 










4 X 4 X 5mm voxel 


21585 


29.4 


-0.54 


63.3 


6 X 6 X 6mm voxel 


7727 


20.1 


-0.46 


37.4 


12 X 12 X 12mm voxel 


980 


9.97 


0.42 


7.01 


Region-based network 


98 


4.49 


1.93 


1.54 



can be erroneously fitted to other distributions fv^ . 

Interestingly, the characteristics of the exponentially 
truncated power law distribution are different between 
the region-based and voxel-based networks. In particu- 
lar, in a voxel-based network, the degree distribution has 
a form P{k) ex ex.-p{—k/kc) with < 0, whereas in a 
region-based network, /3 > 0. In a voxel-based network, 
k^ decreases as k increases, and the decay is further ac- 
celerated by the exponential truncation exp(— /c/Zcc). On 
the other hand, in a region-based network, k^ increases 
as k increases, but this increase is attenuated and even- 
tually overcome for k ^ kchy the exponential truncaton 
exp(— /c/Zcc). If = 0, then the distribution is an expo- 
nential distribution P{k) oc exp(— /c//cc) characterized by 
a single cut-off parameter kc. 

Since the sign of the exponent /3 is the only distinc- 
tion between the degree distributions of the two types 
of networks, it is plausible that the degree distribution 
of a brain network is from a continuum of exponentially 
truncated power law distributions, with P dependent on 
the scale at which the network is formed. To verify 
this, the voxel-based data for one of the subjects were 
down-sampled to larger voxel sizes (6 x 6 x 6mm and 
12 X 12 X 12mm) and the corresponding voxel-based net- 
works were formed at these coarser resolutions. FIG [3] 
shows the degree distributions for voxel-based networks 
at different voxel resolutions, as well as that of the region- 
based network. They all follow exponentially truncated 
power law distributions. Table 1 shows the parameter 
estimates of /3 and kc of P{k) cx k^ exp{—k/kc) at differ- 
ent network scales, as well as the number of nodes and 
the average node degree. As the network resolution in- 
creases from a coarse macroscopic representation of the 
region-based network to a finer mesoscopic representa- 
tion of the voxel-based network, the exponent P decreases 
from positive to negative. The cut-off parameter kc on 
the other hand increases as the number of nodes increases 
in finer representations of the network. Our finding with 
the 12 X 12 X 12mm scale network is consistent with a 
previous study of region-based networks with 1,000 nodes 
demonstrating a clear characteristic of positive /3 in the 
distribution (see Figure SI of [7]). 

Our results above indicate that the scale at which the 
network is formed is important. The higher the resolu- 
tion, the smaller P becomes, resembling more closely a 
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FIG. 3: Node degree distributions 1 — F(k) of voxel-based 
networks of a single subject generated with different voxel 
sizes. The degree distribution of the region-based network for 
the same subject is also shown. 



power law distribution P{k) oc k~^ . This means that the 
network resembles scale-free networks more as its scale 
becomes finer. This may be because connections and 
edges that are represented by a single edge or a single 
node in a coarse network can be represented with mul- 
tiple edges and nodes in a finer network. A single low 
degree node of a coarse network may be represented by 
distinct low degree nodes in a finer network. Similarly, 
the edges at a hub in a coarse network may be resolved 
as collections of edges in a finer network, increasing the 
degree at the hub dramatically. A large number of low 
degree nodes and a very small number of extremely high 
degree hubs are characteristics of scale- free networks [15| . 
It is important to note that p in an exponentially trun- 
cated power law distribution does not become smaller 
than —1 as in a power law distribution of a true scale- 
free network P{k) oc 7 > 0. However, it is possible 
to approximate a power law distribution by modeling the 
distribution with two curves: P{k) oc k~^ for k < kc and 
P{k) oc k-^exp{-k/kc) for k > kc In such a 

parameterization, the tail above the cut-off parameter kc 
is modeled as an exponentially truncated power law dis- 
tribution showing accelerated decay, but the distribution 
of the vast majority of nodes, nodes with degree k < kc^ 
is modeled by a power law distribution. 

It is interesting to note that exponentially truncated 
power law distributions have been observed not only in 
the brain network but also other types of networks, in- 
cluding ecological 0|, social [21, 22], and technological 
0, ^] networks. Such networks may occur due to lim- 
itations in processing capability of nodes or restrictions 
in network growth. Mossa et al. [l9|] generated networks 
with information filtering, or a limitation on nodes' infor- 
mation processing capability, and observed a truncation 
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in the node degree distribution. Similarly, Amaral et al. 
[22] speculated that the cost of additional connections 
to a node functioning near its capacity may limit the 
occurrence of mega-hubs commonly seen in a scale-free 
network. A scale-free network can be formed as a result 
of a process known as preferential attachments [lH, in 
which a new node introduced to a network preferentially 
attaches itself to a high-degree node, resulting a very 
small number of mega hubs. However this preferential 
attachment process may be realistic only when the net- 
work growth is unrestricted. Such a condition is not fea- 
sible for the brain network since there are physiological, 
anatomical, and functional constraints. However, there 
may be some advantages for the brain network not being 
scale- free. A scale- free network is an ideal structure for 
efficiently transmitting information to the entire network, 
but this also means unnecessary or unwanted information 
can be spread throughout the network very easily. Such 
an uncontrollable spread has been observed in epidemi- 
ological networks in which a disease spreads to the en- 
tire network without any threshold for epidemics (2^, [^J] . 
On the other hand, an exponentially truncated power law 
network is resistant against such massive epidemics, with 
a non-zero threshold for an epidemic [IS] . In the context 
of the brain network, a scale-free structure could cause 
extensive synchronization often seen in epilepsy patients, 
while an exponentially truncated network may prevent 
such over synchronization. This may also allow localized 
processing in a certain module of the brain without in- 
volving the entire brain network unless the demand of a 
cognitive process exceeds a certain threshold. 



In summary, degree distributions of macroscopic and 
mesoscopic functional brain networks were examined, 
through constructing region-based and voxel-based net- 
works, respectively, on the same subjects. It was found 
that the degree distributions of both networks follow ex- 
ponentially truncated power law distributions, and not 
power law distributions as previously reported on voxel- 
based networks. The discrepancy in the literature is 
likely due to the fact that prior studies of voxel-based 
networks plotted the degree distributions P{k) whereas 
all region-based network analyses plotted the cumula- 
tive distribution 1 — F{k). Due to sensitivity to noise 
in the extreme tail of the degree distribution, the cumu- 
lative distribution plots are preferred. It was also found 
that the distributions of region-based and voxel-based 
networks belong to the same continuum of exponentially 
truncated networks, and finer scales of the network result 
in networks more similar to scale- free networks. The re- 
sults from this study suggest the importance of modeling 
the brain network at resolutions as fine as data would 
allow. 
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